Kondo effect in complex mesoscopic structures 
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We study the Kondo effect of a quantum dot placed in a complex mesoscopic structure. Assuming 
that electronic interactions are taking place solely on the dot, and focusing on the infinite Hubbard 
interaction limit, we use a decoupling scheme to obtain an explicit analytic approximate expression 
for the dot Green function, which fulfills certain Fermi-liquid relations at zero temperature. The 
details of the complex structure enter into this expression only via the self-energy for the non- 
interacting case. The effectiveness of the expression is demonstrated for the single impurity Anderson 
model and for the T-shaped network. 



I. INTRODUCTION AND OUTLINE 

The single-impurity Anderson model (SIAM) has been 
a paradigm of a strongly correlated electron system. 1 
This seemingly simple model gives rise to dynamical 
screening of the local spin by the electrons in the Fermi 
sea, leading to a crossover from a weak coupling system 
at high temperatures (T) to a strongly coupled one at 
low T, with the relevant temperature scale given by the 
Kondo temperature, Tk- The recent observation of the 
Kondo effect in quantum dots (QD's), whose parameters 
can be tuned continuously, and which allow probing of 
various properties, 2 has yielded strong theoretical efforts 
in this direction. 3 Recent experiments on a QD embedded 
on one branch of the Aharonov-Bohm interferometer 4 
created additional interest in Kondo effects in complex 
networks. 

While the high-T behavior of the SIAM can be ade- 
quately described by perturbation theory or poor-man 
scaling, 5 and the low-T behavior is described by Fermi 
liquid theory, 6 there is no simple theory that describes 
the model's dynamical properties correctly in the whole 
temperature range, including both the high- and the low- 
T limits. This crossover has been described by the 
computationally demanding numerical renormalization 
group (NRG) approach. 7 Other methods to describe the 
crossover are the "conserving T-matrix approximation" 8 
(which overestimates the unitarity sum-rule), or the more 
limited "non-crossing approximation" , 9 and to some ex- 
tent, quantum Monte Carlo calculations. 10 Some of 
these methods have also been applied to the QD in 
the Aharonov-Bohm interferometer. 11 However, none of 
these methods has the flexibility to follow analytically 
the effects of various network parameters (e.g. a mag- 
netic flux) on experimentally measurable quantities, in 
particular dynamical ones. 

In this paper we discuss the SIAM for a QD which 
is embedded in a general complex network. Many of 



the interesting physical properties of the system can be 
expressed explicitly in terms of the single electron (re- 
tarded) Green function (GF) on the dot, Gdd<r(u), for 
electrons with spin component a and energy u> (measured 
relative to the Fermi energy). Here, 



G ABa {w)=<^ A a \B\ > 



(1) 



(For the dot GF, A = B = d a , where d a destroys an 
electron with spin a on the dot). One example concerns 
the local density of states (LDOS) on the dot, 



Pda{u) = -3[G ddcr (w)]/7r. 



(2) 



Another example is the conductance Q between two leads 
which are connected to sites on the network. For a gen- 
eral mesoscopic structure, this (linear response) conduc- 
tance has the form 12 



2e^ 

h 



/<K-tD», 



(3) 



where f(ui) is the Fermi distribution function (we set the 
Fermi energy at ui = 0), 



/(«>) = 



1 



oiu/kT _|_ I ' 



(4) 



and G(uS) is often given in terms of the retarded Gdda-, 
the advanced G* dda and non-interacting parameters char- 
acterizing the network. For the simple SIAM (dot and 
leads), Q{oj) is proportional to pd CT (o;). 13 ' 14 

Here we use the equations of motion method to derive 
a simple analytic approximate expression for Gdda, for a 
QD on a general network. Our approximate formula ex- 
hibits the correct behavior both at high and at very low 
temperatures. For simplicity, we consider a QD with a 



single level (of energy e<j)- Electron-electron interactions 
are assumed to exist only on the QD, and we take their 
energy U to be infinite. In our scheme, one first solves for 
the GF on the dot in the absence of the electron- electron 
interactions, i.e. U = 0. This (spin-independent) GF can 
be written in the form 



GSdM 



1 



+ in-e d - S (w)' 



w 



(5) 



with the energy dependent (complex) self-energy (SE) 

Z Q (u) = 6e d (u>)-iA(u>). (6) 

This is easily done exactly, for any finite network; it only 
involves the solution of a finite set of linear equations. 
Our approximate explicit expression for Gdda{u) is then 
given in a single equation [see Eq. (53) below]. Some- 
what surprisingly, this equation depends on the network 
parameters only via Eq(w). Although approximate, this 
expression allows for detailed systematic investigations of 
the Kondo effect as function of the system parameters, 
on a variety of complex networks. Even if sometimes 
only qualitatively correct, such systematic studies help 
to investigate new physical phenomena on a broad vari- 
ety of mesoscopic systems. Furthermore, the dynamical 
mean field theory 15 (which was developed to address the 
physics of the periodic Anderson and Hubbard models) 
iterates the local GF of the SIAM, which is calculated at 
each stage in terms of the effective SE created by the rest 
of the lattice. Since our GF is easily calculated in terms 
of E , it is an ideal candidate for such calculations. We 
are not aware of alternative simple analytic expressions 
which obey the necessary requirements at both high and 
low T. 

In Sec. II we follow Rcfs. 16-19, and derive the 
equation-of- motion (EOM) for the dot GF. Extending 
these references, this EOM is derived here for a gener- 
alized case, in which the dot "sits" within an arbitrary 
complex network. The EOM for the dot GF involves 
higher-order GF's (consisting of more operators), whose 
EOM's introduce in turn more GF's. One then termi- 
nates this hierarchy by decoupling out averages of oper- 
ators. The latter are then found using the fluctuation- 
dissipation relationship, from the relevant GF's. In this 
manner, the treatment becomes self-consistent. The suc- 
cessful decoupling must keep as much of the electronic 
correlations as possible. 17 For example, a widely used ear- 
lier approximation 14 neglected some dot-lead correlations 
and therefore gave reasonable results only at T > Tr-. 
For this reason, Gcrland et al. had to combine the 
EOM method at high T with the NRG at low T. Af- 
ter we correct these earlier calculations, and include all 
the necessary correlations, we obtain good results for all 
T. An earlier application of the EOM method solved the 
integral equations numerically, on the simplest interfer- 
ometer geometry. 21 

The above decoupling scheme produces an integral 
equation for Gddo-(w), which generalizes that found in 



Refs. 18 and 22. A new approximate analytical solu- 
tion for this equation is found in Sec. III. The low 
temperature limit, the Fermi liquid conditions, and the 
Kondo behavior are discussed in Sec. IV. Unlike earlier 
papers (spanning over four decades), our solution gives 
good qualitative results for the whole parameter range: It 
has the correct high-T behavior and fulfills several low-T 
Fermi liquid relations. 

For the simplest SIAM, the "network" is represented 
by a band with a density of states N(ui). In the broad 
band limit, the self-energy on the QD can be approxi- 
mated by its value at the Fermi energy, Yiq{uS) w —iAt = 
—iwV 2 N(0), and A& represents the width of the impurity 
state in the absence of interactions (V represents the ap- 
proximately energy- independent dot-band coupling). In 
this case, our equations reduce to those discussed in the 
earlier literature. We emphasize again that even the so- 
lution for this simple case is new and non-trivial, since 
earlier work cither did not have a close analytical formula 
or missed some correlations, leading to wrong results at 
lowT. 

As stated, our result for Gdda{u) depends on the de- 
tails of the general network only via the non-interacting 
self-energy, Eo(w). To demonstrate the power of our re- 
sult, we present here a few simple examples. Section V 
contains explicit results for the simple case of the SIAM, 
when the QD is coupled to two one-dimensional leads. 
For this case, we easily calculate various features of the 
Kondo behavior, including the peak in pdda at the Fermi 
energy, the plateau in the conductance in the so-called 
unitary limit and the plateau in the "transmission phase" 
at tt/2. 

A second example, of the so called "T" network, is 
solved in Sec. VI. In this case, one has interference be- 
tween the wave function on the dot and that on the inter- 
section point. Our simple approximate formula exhibits 
the Fano vanishing of the transmission and the associated 
so-called anti-Kondo effect, as seen experimentally 23 
Here we also extend earlier theoretical work, which used 
complicated techniques. 24 ' 25 Finally, Sec. VII contains 
our summary. 



II. EQUATIONS OF MOTION 
A. The Hamiltonian 

Our Hamiltonian, 

n = ^dot + Wnct + ^dot-not; 

contains the dot part 

Ti-d = £d ^2 n da + Un d \n di , (8) 



(7) 



with a single energy level q and an electron-electron in- 
teraction U (nda — o%d a ). All the other parts of H can 
be arbitrary, but non-interacting: 



H u 



E e « a * 



E J ' 



7im Q* n Q-Clrna 



describes the network without the dot, and 



Wdo 



2_^(Jdndla r , 



h.c. 



(9) 



(10) 



describes the coupling of the dot to the network. Here, 
a nt7 creates an electron with spin a and energy e„ on 
the site n of the network. The coefficients J nm and Jdn 
represent tight-binding hopping matrix elements. In the 
absence of a magnetic flux, these coefficients can be cho- 
sen to be real. The flux then turns them complex, with 
phases which relate to the Aharonov-Bohm effect and 
with J nm = J^. 11 - 21 - 26 

The main part of this paper concerns the above gen- 
eral Hamiltonian, which assumes no specific details on 
the structure of the network. However, for some specific 
applications, it is convenient to attach a finite mesoscopic 
network to several semi-infinite leads, which connect to 
electron reservoirs. We thus divide W nc t into two parts, 



H, 



H n 



Tii 



cads; 



(11) 



where 7i mcS o has the same form as H nc t, except that the 
sum is restricted to sites which belong to the finite meso- 
scopic part of the net (excluding the dot and the leads) , 
and 

Hlcads = J2i HL - E (^.'C^' + h - C )\ > 
L n L ,r 

H L = -J2JreC& )1 c'£>+h.C.. (12) 

rs 

Here ul are the indices of the points on the mesoscopic 
structure which are connected to lead number L, while 
c ra destroys an electron of spin a on the site r of the 
lead. 

Usually, the lead L has a continuous spectrum, with 
eigenstates |fc) and eigenenergies €k, in a band of width 
2D. Equation (12) can then be written as 

w lcads = E [Hl - E ( y «- (*)«U4? + h - c )} > 

L til ,k 

<kCka" C ka> (I 3 ) 



n L = Y l -" {L) ^ {L) 



and also 

't-dot — net / j \Jdn@'(jQ"n<j 

na 

-EE^W^ + kc), (W) 



L k(T 



with obvious definitions. When the network consists only 
of the dot and the leads (i.e. H m0 so = 0) then this has 
exactly the form of the standard SIAM model, discussed 
broadly in the literature. If all the leads have similar 



bands, then the electrons in the whole system also have 
a continuous spectrum in the range —D<lu<D. All 
integrals over u> below will be thus over this energy band, 
and it is due to resonances in such integrals that we need 
to calculate the retarded GF at u + in. 



B. Derivation of Eo(w) 

From now on we return to the general Hamiltonian, 
Eqs. (7-10). The EOM for the retarded GF on the QD 
reads 

(w + in - e d )G d da = 1 - E J dnG nda + UT a , (15) 

n 

where 

r„=«n,i_A;4»- (16) 

From now on, we shall use the shorthand u for u> + if], 
except where it matters. Specifically, 



\U) — e n )G n da — —Jdn^dda ~ / J Jnm.G 



rada • 



(17) 



We next define the matrix 

M(u)) nm = (n\u> - Hnet\m) = Jnm + 5 nm {uJ - £„), (18) 

which represents the system without the dot. In the fol- 
lowing, we shall need the inverse matrix, 



T{u) = [M(w)\- 



(19) 



which is the Green function for the system without the 
dot. With this definition one finds 



where 



G m da — U m {UJ)&ddm 



,(w) = -E-^MrrmJn 



(20) 



(21) 



We can now calculate the second term on the RHS of Eq. 
(15), and show that 



/ J JdnGnda — SqG ddcr , 



(22) 



with 



^o(w) = — 2_^, JdmUmiuj) = ^ Jdm,J~{u)mnJnd- (23) 



At U = 0, Eq. (15) thus reduces to Eq. (5). Clearly, E 
can be calculated for the non-interacting case. 

All the other GF's can similarly be expressed in terms 
of Gdda- For example, 



Gdmcr — u m ( w )Gddo 



(24) 



where lu* = uj — irj. Another way to represent u^(oj*) is 
to start from u m (u)), and then take the complex conju- 
gate only of the hopping coefficients J mn (namely reverse 
all their phases), without touching ui + irj. Similarly, 



Glma — J-(uj)i m + Ue(u})Gdmcr- 



(25) 



Since M. is an infinite matrix, its inversion may not be 
trivial. Formally, we can denote the eigenenergies and 
eigenstates of 7i ne t by e a and \a), and then write 






1 



IT] 



-(a\n). 



(26) 



However, in the examples it is useful first to eliminate 
the parts which involve the sites on the leads. Examples 
for this procedure are presented in Appendix A. 



C. EOM's for higher-order GF's 

The EOM for r CT (with only on-site Hubbard interac- 
tion on the dot) does not involve GF's of more arguments: 

(w - e d - U)T a = (n d ^) + ]T J* dn T^ 

n 
n 

with three new GF's, 

r ni = < nd-<Ta na ;dl >, 
TV* = « ($_ a a n - a d a ; 4 », 

ri 3 J = «4-^-^;4»- (28) 

The inhomogeneous term in Eq. (27), {nd~ a ), represents 
the average number of electrons with spin a on the dot. 
This number needs to be determined sclf-consistently, via 
the fluctuation-dissipation theorem: 



(nda) = / duf(ui)pda(u>)- 



(29) 



We shall return to this condition below. 



W 



We now write the EOM's for the new GF's, Tna, and 
use several approximations for their solution. Firstly, it 

(3) 

is easy to see from the EOM for T n< j that this GF is of 
order 1/U. Therefore, it contributes to r^ only at order 
1/U 2 , and can be ignored for U — ► oo [we need T a only 
to order 1/U, see Eq. (15)]. Secondly, we introduce a 
decoupling scheme for the new GF's that appear in the 
EOM's of there's, 

^ ^n<ri @"wi(T2 ^c"3 5 a ■^ > — y^na-i ^mai )Q<Ji<J2 < ^ ^o^i^o -^ 

- («Li^i)^i» 3 < a m<J2 ;dl > . (30) 



/I / 



We omit thermal averages of the form {a na d a ), which 
include two destruction (or creation) operators, and 
(aJj CT a mCT ') with a ^ a', relevant only for states with a 
net magnetic moment. The latter assumption means that 
we restrict the discussion only to symmetric states, with 
Gddi = Gddi- 

Using these approximations, the EOM for r„ „ becomes 



(w - e n )T^ = -Jdn T ° - 



E 



j r (1) 



(31) 



The derivation of this equation also required an ad- 
ditional term, -G n dcr[E m J dm (d' ! _ <T a m - a ) - c.c.}. How- 
ever, this term vanishes. Here and below, we calculate 
equilibrium thermal averages by using the fluctuation- 
dissipation theorem, 

{At B) = / (" S /M ( GBA G * ab ) ■ (32) 

Equation (31) is practically the same as Eq. (17), 
which was used to derive So in terms of the J nm 's. Ap- 
plying the same algebra to Eq. (31) then yields the ana- 
log to Eq. (22), 



2 , Jdn^nl = Sor CT . 



(33) 



^(2) 



Using the same approximations, the EOM of r„<j is 
(t i - e \r( 2 ) — - t* r _ V^ t r( 2 ) 

yu/ t n )L na — '->dn L a / , • J nm i - mG 

rn 

+ {X n )G dd a + (dl a a n - a ) (1 + Z Gddv) , (34) 

in which 

(X n ) = E Jdm( a ln-a a n-a) ■ (35) 

m 

Although the first two terms on the RHS of Eq. (34) are 
similar to those in Eq. (31), we now need to calculate the 
other two terms. Since these terms require thermal aver- 
ages, which we express using the fluctuation-dissipation 
Eq. (32), they end up with integrals which involve Gdd-o, 
leading finally to our integral equation for the dot GF. 
The details of these calculations are presented in Ap- 
pendix B. 

Inserting the results for the necessary thermal averages 
into the RHS of Eq. (34), the results are used in Eq. (B5) 

for the combination — ^2 n Jdn^na which is needed in the 
RHS of Eq. (27) . Adding also the corresponding combi- 
nation for r^ 1 -* [Eq. (33)], yields our final result for T a . 
It remains to insert it back into the EOM (15) for the dot 
Green function Gdda- This yields an integral equation for 
this function, 



du>'\ ,/V) 



2iri) ui — uj' + irj 



(1 



^oGdd- 



r )(£ -£o)-(u/-u/*) 



(1 + EqG^) I [-P, 



2iriJ lu 



/(«') 



ITj 



G^Eq - E ) - (u/ W*) 



r/^O 4 



(36) 



Here and below, G dda and Eo are understood to be functions of uu — > a; + irj, and the primes denote a dependence on 
u/ — ► u/ + i?j, with u/* — > u/ — ir/. Two comments are in place here. First, for the simplest SIAM, when Eo(w) ~ — JA5 
is independent of to, we have Eo(u>*) « iAj, so that only the parts with uu'* survive. Furthermore, in this case one 
can factorize Aj, out of the integrals. This reproduces the integral equation of Lacroix, 18 



(w - e d + iA b )G d da = 1 - (nd-a) - 2iA(, 



(l-iA b G Ma ) l(-%L 






27T«/ LU 



ir\ 



\G'a 



dd—a\ 



(37) 



However, even for this simple case there has not been an 
analytic solution that covers the whole parameter range. 
Second, discarding 14 correlations between the dot and 
other sites on the net, e.g. (dj.o no -), amounts to neglect- 
ing the terms containing G' dd in Eq. (36). This ends 
up with a breakdown of the Fermi liquid conditions at 
T = 0, and with a bad approximation for T < Tk- 



III. APPROXIMATE SOLUTION OF THE 
INTEGRAL EQUATION 

We now restrict the discussion to "non-magnetic" 
states, and replace Gddi = Gdd{ = Gdd- Equation (36) 
can then be written as 

(w - e d - >Jo)Gdd = Sn + AG dd - B(l - ^AG (W ), (38) 



where 5n = 1 — (n^) = 1 — ( n di) = 1 — 
and B arc functions of u>, given by 


(nd)/2 ar 


A 


/■/ duA /(«*) 




7 V 27ri J uu — uu' + irj 


X 


s o (i + G dd( s o - S ) J - (u/ + ir? -> u/ - ir?) 


B( 


/■/ <fcA /(u/) 




J V 27rz / uj — uj' + zr/ 


X 


GI w (EJ,-Eo)-(u/ + m?-c/-m?) 





(39) 



For reasons that will be explained later, it is now con- 
venient to replace /(u/) inside the integrals by [/(a/) — 
1/2] + 1/2. The part related to the constant 1/2 can then 
be calculated using the Kramers-Kronig relations, 



du>' 



V 



1 



2iri I uj — uu' + irj 
and thus we find 



(*V)-*V)*)=*V)> (40) 



^(W) = -E W+ (-7T-) '■ 

2 J V 2m / uu — uu' + rq 



S (l 



(l + G^(E - E )) - (u/ + 27? - u/ - «,) 



B(w) = 



duA /(q/)-i/2 

27ri/ w — lu' + irj 



G^(E -Eo)-(u/ + ir ? ^u/-ir ? ) 



(41) 



We next follow Lacroix, 18 and assume that the integrals 
are dominated by the region uu' ~ w, namely that 



/ duu'F(uj, uj 



,J(u>')- 1/2 



lu — lu' + irj 
and thus 

B(w) « A(u,)G^(u;)X(a;), 
where A(w) = — 3E (w) and 22 



F(w,w)X(u;), 



(42) 



(43) 



X(w) 



13 <L/ ,/V) - 1/2 

2 ln (2.) 2 "*<; 



27ri ; 



(44) 



(\I/ is the Digamma function). A similar procedure yields 



A(lu) « A(w) [-- + X(w)(l + »A(w)G*(w))] . (45) 

The reason for using the above transformation is that 
when one analyzes the equivalent equations for finite U, 
then one requires particle hole symmetry, namely the re- 
lation 



G M (u) = -G dd (-cu), 



(46) 



where the 'tilde' denotes particle-hole transformed quan- 
tities. In that case, this relation is equivalent to the re- 
lations 



B(lu) = -B*(-«), 

A(cu) = A*(-oj)-^* (-uj). 



(47) 



These relations should hold also for any approximate 
solution. 16 It is easy to check that these relations are 
obeyed by our approximate expressions (43) and (45), 
but would not hold if we made the approximation (42) on 
the original equations, before shifting f(uu') by 1/2. 16 ' 18 



Taking Sn as a parameter, and defining the scaled vari- 
able 



z(u>) 



u - td - Se d (ui) 



2A(w) 
the solution of Eq. (38) can be written as 

G dd {uj) = g{u>) (Sn + iQ(u>)/X*(v 

with 



1 



9 



u> - e d - 5e d + i'iA/2 2A(z + 3i/4) : 



Q = S-[S 2 + \X\ 2 q5n-6n 2 )} 1 / 2 , 

S = z 2 + -^ - zMX + (Sn - h$!X, 
16 A' 



(48) 

(49) 

(50) 

(51) 
(52) 



and we have omitted the explicit dependence on u for 
brevity. For some purposes it is more convenient to write 
Eq. (49) as 



G dd {uj) 



1 



1 - C(w) u;-(e d + SE d (u)) + t|A(w) ' 



(53) 



where the real functions C(tv) and SE d (uj), which result 
from the strong interactions on the dot, depend on T and 
on Eq(u') and are given by 



8E d {uu) = Se d - 2A(KX)C/(1 - C), 

CM = 1 = 2(z 2 + 9/16) 

2A 2 g| ff (z + (5nX*/Q)| 2 Q|* + <JnX*/Q| s 



(54) 



Equation (53) forms our general result. We emphasize 
again that all we need to know is the non-interacting 
self-energy, S (w). 

As seen explicitly from the above equations, our ana- 
lytic expression for G dd depends on the parameter Sn — 
1 — (n<j)/2, which needs to be determined self-consistcntly, 
via Eq. (29). We discuss this condition in a specific ex- 
ample below. 



IV. THE FERMI LIQUID RELATIONS AND THE 
KONDO TEMPERATURE 



At T = 0, one easily sees that 

SRX(w) w In |D/w|/7r, 



(55) 



diverging at the Fermi energy oj — > 0. At T = ui = we 
thus have C — > and C5RA — ► u(z), with 



2(z 2 + 9/16)[z + Vz 2 + fa(1.5 - fa)] 

IMZ) = ; , (56) 

[z+^z 2 +Sn(1.5-Sn)] 2 +5n 2 



and z = z(u> = 0) = -[e s + 5e d (0)]/[2A(0)}. Therefore 

1 



G dd (0) 



A(0) ((4/3) [«-«(«)]+») 



(57) 



For a very deep level on the dot, z»l, u(z) approaches 
z, and G dd approaches — i/A(0), ending up with high 
plateaus in the LDOS and in the phase of the complex 
G dd {uj) sticking to tt/2. All of these characteristics, which 
are hallmarks of the Kondo behavior, 20 are observed in 
the examples below. As e d moves up from large negative 
values, Eq. (57) shows a gradual crossover away from 
this "unitary" limit, and G dd becomes mainly real and 
small as z — > — oo. 

Eq. (57) also ensures that the GF obeys the Fermi- 
liquid relations 27 at T — 0: (i) The imaginary part of 
the SE on the Fermi level is the same as in the ab- 
sence of interactions, i.e. — A(0), confirming the uni- 
tarity limit; 28 ' 13 (ii) For e d <C -A, the LDOS there 
[i.e., p d {u) = 0) = pd<r(0)] approaches the finite value 
l/[7rA(0)]. Wc emphasize that these results apply for 
any complex network; the only input is £o(w). 

The Kondo effect concerns the behavior of the density 
of states near the Fermi energy at low temperature. In 
that regime, G dd is dominated by its imaginary part, and 
thus 



W «1/[ttA( W )(1-C)]. 



(58) 



For diX(uj) > z(lu) > 1, rne has 5 « z(z - MX) < 0, 
Q w 25 and the leading w-dependence of C is 



C 



z-StX 



(59) 



We thus end up with a logarithmic cusp in p d at u> = 0. 
The related narrow peak reaches one half of its peak value 
(equal to l/[7rA(0)]) when C w — 1. One definition of the 
Kondo energy Tk is to identify Tk with the half-width 
of this peak. Solving C(uj — ±Tk) = — 1 and using Eq. 
(59) thus yields MX = 2z, i.e. 



T K (ed) = Dex.p[-air\e d + 5e d {0)\/A(0)}, 



(60) 



with a = 1. This result, which agrees with that 
of Lacroix, 18 is qualitatively similar but quantitatively 
different from the presumably exact Tk as given by 
Haldane, 29 which has a = 1/2. However, Tk only rep- 
resents some crossover energy scale, and we expect the 
above solution to follow the qualitative variations of Tk 
with the system parameters, which are contained in Sed 
and in A. It should be emphasized that Tk depends 
on both A and Sed (some authors ignore the part com- 
ing from the real part of the SE, <5e<j). For example, 
for the Aharonov-Bohm interferometer Tk oscillates with 
the flux. 30 

At the Fermi energy, ui = 0, and for fixed ed < 0, one 
expects pd{0) to decrease with increasing T. In this limit, 
X(lu = 0) = ln|AD/T|/7r, with A = e -*( 1 / 2 )/(27r) « 



1.13387. 
relation 



For \ed\ 3> A(0) we then find the asymptotic 



MO) 



1 



MX 



ttA(0)(1-C) 
1 /- . 7r(e<i 



ttA(0)KA 



ttA(O) V 2A(0)ln(AD/T)/' 



(61) 



This density of states decreases to one half of its maxi- 
mum T = value at a temperature ATk , with Tk given 
in Eq. (60). Although the result (61) is qualitatively 
correct, and is also consistent with our Tr-, its explicit 
low-temperature dependence disagrees with exact expec- 
tations: we obtain a logarithmic behavior, whereas Fermi 
liquid theory predicts a T 2 dependence at low T. 6, This 
is an artifact of our approximation. However, at high T 
we do recover the usual logarithmic variation. 

In the opposite limit, of \z\ ~3> KA ~ 1, one finds that 
\Q/X\ <C 1, and Eq. (49) implies that Gdd{u) approaches 
Sn g(u>). The appearance of Sn in the numerator, in 
place of 1, results from the fact that some weight of the 
spectral function is pushed to infinite [7. 14 Another in- 
teresting point concerns the factor 3/2 in Eqs. (50) and 
(53). This factor may be explained heuristically 31 by the 
fact that while both spin directions are accessible for tun- 
nelling into the dot, only a single one can tunnel out of 
it. 



V. EXAMPLE 1: THE SIMPLE SIAM 

In the following two examples, we assume simple semi- 
infinite one-dimensional leads, with identical nearest- 
neighbor hopping matrix elements, —J, and with a lat- 
tice constant a. The eigenenergies of each TCl are 
therefore tk = — 2 J cos ka, with eigenfunctions (n\k) = 
sinnkay / 2/£l (tt — > oo is the length of the lead), and 
D = 2J. We also assume that the leads L and R are 
attached to the mesoscopic network only at one site, 
with hopping matrix elements Ji and J r . For the simple 
SIAM, when the two leads are directly connected to the 
dot, we show in Appendix A that the leads generate a 
self-energy 

Eo(w) = -e ik|a (l M 2 + \Jr\ 2 )/J = -e^ a A b (62) 

on the dot, with q determined by u> = — 2 J cos qa and 
with 



A 6 = r^ + r r . 



Here we have used the notation 



|J,, r | 2 /J = ^(0)|J A ,, 



(63) 



(64) 



with Af(u>) being the density of the band states. The 
width of the non-interacting resonance at e<j is thus 
equal to Af,smqa. In the figures presented below we 
use the symmetric case, Tg = T r = A&/2. Substituting 



S (w) = Se d (uj) - iA(u) into Eqs. (49)-(52) then yields 
Gdd- 

We start by discussing the self-consistency condition 
for (rid). Figure 1 shows an example of our self-consistent 
solutions for (rid), at T = 0: we start from an initial guess 
for Sn, then calculate (rid) from Eq. (29), and iterate; the 
procedure converges after a few iterations. Also shown 
is n, which represents the total change in the electron 
occupation in the system due to the dot, as determined 
from the Friedel phase, defined via 13 ' 28 ' 27 



tan[-(l 



n)}=^[Gdd(0)}/^[G d d(0)}- 



(65) 



For ed <C — A(0) = — Af,, n approaches 1, while (rid) re- 
mains slightly smaller. This small difference could reflect 
an additional occupation of other sites in the leads. It 
could also result from the inaccuracy of our approxima- 
tion for Gddis^), which becomes worse as u) moves away 
from the Fermi energy; the integral in Eq. (29) contains 
contributions from all u). However, this small difference 
has only a small effect on the other calculations presented 
below. 




FIG. 1. The self-consistent average occupation on the dot, 
(rid) (thick line), and the Friedel added occupation h at T = 0, 
as function of td, for the simplest SIAM with Ab/J — 0.1. Re- 
sults remain almost unchanged up to T ~ Aj,. 

Although it is easy to solve for (rid) for each set of pa- 
rameters and use the resulting 5n for other calculations, 
the qualitative results are only weakly affected if one uses 
an arbitrary smooth variation of (rid) from to 1 as td 
varies from +oo to — oo. In the calculations below, we 
used such a smooth variation for Sn, with a width of 
order A&. 

An alternative approximate expression for Sn follows 
from Eq. (49). It turns out that for most of the inte- 
gration range in Eq. (29), Gdd(w) is dominated by the 
first term there. Neglecting the second term, one has 
1 — Sn = (rido-) ~ Sn n$, with no = — / <ko'^[g(uj)]/TT de- 
pending only on the non-interacting parameters. Thus, 
Sn w 1/(1 + no). For A;, between 0.1 and 0.01 and for 
£d < Ah, we find this estimate to be within ~ 10% from 
the full self-consistent value. 

Figure 2 presents the LDOS [Eq. (2)] versus u), for the 
simple SIAM, with the parameters as indicated. Plots 
at lower T are indistinguishable from the one shown at 



T = .01J. Our approximation reproduces the Kondo 
peak at low T, in addition to the much broader peak at 
td- Figure 3 shows pd at the Fermi energy (oj — 0) versus 
Ed- Note the very slow (logarithmic) increase of pda with 
decreasing T, in agreement with Eq. (61). 
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FIG. 2. The LDOS on the dot, for the simple SIAM, with 
ed = —.3, Ab = .1 and T = 0, .01, .05 (all energies are in 
units of J). The Kondo peak values near u) — 0, shown with 
more detail in the lower panel, decrease with increasing T. 




FIG. 3. The LDOS at the Fermi energy, pd{0), versus td, 
for T = 0, 10~ 30 , 10~ 10 , 10 -3 and A b = .1; p d {0) increases 
with decreasing T. 

We now turn to the conductance, Q. For the simple 
SIAM, one has 13 ' 14 



Q(u>) = 



4T f T r 



■SGdd{u)- 



Equation (53) implies that 

SGdd = -(l-C)A b \Gdd\ 2 . 



(66) 



(67) 



For T = 0, when C — ► 0, we also have Q(u>) = 
AYzY r \Gdd\ 2 - At finite temperatures, these two expres- 
sions for the conductance exhibit the same qualitative 



behavior as pd (Fig. 3). The quantitative difference be- 
tween the two expressions, which increases with T and 
with large negative e<j, represents a breakdown of unitar- 
ity which may be an artifact of our approximation. 

At low T, (—df/dco) is practically a delta function, 
and Eq. (3) yields Q « (2e 2 /h)G(0), exhibiting the qual- 
itative behavior shown in Fig. 3. At higher T, the peak 
in Q slightly below td — becomes lower and broader 
than that of pd(to = 0). At fixed td < 0, increasing T 
results with a decreasing Q, but with an interesting su- 
perimposed peak at T = 0{\ed\) (when the peak at e^ 
starts to contribute). 

Some experiments 4 place the QD on one branch of 
an open Aharonov-Bohm interferometer, and attempt 
to extract the transmission phase for scattering through 
the dot. The transmission phase is usually related to 
the phase of Gdd- 20 Without getting into the question 
of what is really measured in the interferometer, 32 it is 
still of interest to study the latter phase. Figure 4 thus 
shows the "transmission", represented by A 2 |G(w| 2 , and 
the "transmission phase", represented by [see Eq. (65)] 
a = 7r(l — n/2) as function of td for several tempera- 
tures. Interestingly, at high T this phase simply grows 
smoothly from zero to 7r through the resonance, similarly 
to the non-interacting case or to the Coulomb blockade 
case. However, as T decreases, the peak in the "transmis- 
sion" broadens towards negative e^, eventually reaching a 
plateau for T = (see also Fig. 3). At the same time, the 
phase develops an intermediate plateau at ir/2. 20 This 
plateau begins at an energy td which is roughly given by 
T ~ Tjc(td), as defined in Eq. (60). Studying the energy 
where this phase grows from zero to tt/2 thus suggests 
another way to define the crossover temperature Tr-. 
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FIG. 4. "Transmission" and "transmission phase" (in 
units of 7r) for the simple SIAM, for Aj, = 0.1 and for 

r = o, lo -10 , io -3 . 



VI. EXAMPLE 2: THE QD ON A "T" NETWORK 

Our equations also become very simple for the "T" 
network, when the dot sits on a side branch. Such a 
network has recently attracted both theoretical 24 ' 25 and 
experimental 23 interest, as the simplest realization of the 
Fano-Kondo effect. The mesoscopic network now consists 
of a single site "0" . This site has one bond connected to 
the dot, with hopping J x , and two bonds connected to 
two leads, with hopping Jg and J r . In this case, the ma- 
trix Fu in Appendix A is of order lxl, and Eqs. (A3) 
and (A8) yield T m = l/[w - e + e l ^ a (\J e \ 2 + \J r \ 2 )/J}- 
Using A;, from Eq. (63), we have 



S (w) = 



'' X 



u - e + e i l'?l a A h 



= 6e d {u) - iA(u>). (68) 



2.1 



In this case, Q of Eq. (3) is given by 

£M = - * SGoo(^). 



r r 



(69) 



The calculation of Goo 
and (25), 

Go 



follows directly from Eqs. (24) 



1 + S G 



dd 



LU - € Q + iA b 



(70) 



Therefore, SGoo is a linear combination of tflGdd and of 
3Gdd, with coefficients which depend only on the non- 
interacting parameters. 

These equations reproduce those found e.g. in Refs. 24 
and 25. However, at this point those authors use com- 
plicated numerical schemes to obtain G^, which do not 
allow for systematic studies of the dependence on the 
various parameters. In contrast, we can easily obtain the 
approximate Gdd analytically, with the same ease as for 
the previous example. All we need to do is to substitute 
Eq. (68) into Eqs. (49)-(52). 

For the non-interacting case, Gdd is simply given by 
Eq. (5). Figure 5 then compares the results for the con- 
ductance (in units of 2e 2 //i), with and without interac- 
tions. Without interactions (thin lines), for eo = the 
upper panel in Fig. 5 shows a symmetric Fano vanish- 
ing of the conductance at e^ = 0. As e becomes more 
negative (lower panel there), the conductance assumes a 
typical non-symmetric Fano shape, with smaller conduc- 
tances at large \ed\- The interactions (thick lines) have a 
negligible effect at ed ^> 0, when the dot is not occupied. 
However, at e^ < and T = (thickest lines) the inter- 
actions cause strong changes. For eo = and large neg- 
ative ed, the occupation on the dot causes the so-called 
anti-Kondo effect, 24 where the formation of the Kondo 
singlet causes destructive interference which yields zero 
conductance. For eo < 0, the Fano zero is shifted to lower 
e<£, the peak disappears and the asymptotic conductance 
increases. As T increases (medium thick lines), the be- 
havior at negative €d gradually returns to that of the 
non-interacting case, and the Fano zero disappears: the 



conductance is always non-zero. Again, all of these phe- 
nomena are qualitatively similar to those found in Refs. 
23, 24 and 25. 
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FIG. 5. Conductance of the "T" network, with 
J x = .2, At = .1 (all in units of J). Thin lines: without inter- 
actions. Thick lines: with infinite interactions at T — (in- 
termediate thickness lines: T = 10~ ). Upper panel: eo = 0. 
Lower panel: eo = — .1. 

A slightly more complex T-network is obtained by plac- 
ing another non-interacting site, "1" , between the site 
and the dot, with energy t\ and with hopping elements 
J x to both site and the dot. Our mesoscopic system 
now contains two sites, and 1, and the methods of Ap- 
pendix A yield 



S n = 



Go 



J 2 (cu-e + e^ a A b ) 



D 



to 



ei 



D 



D 2 



:G 



(Id; 



with 



£> = («- ei)(« - e + e^ a A b ) - J 2 . 



(71) 



(72) 



(73) 



Figure 6 shows a few examples on how the interplay 
between eo and e\ can change the dependence of the con- 
ductance on td for the fully interacting case at T = 0. 
For ei = eo = 0, the graph looks exactly like that for the 
simple SIAM, discussed in Sec. V. Clearly, the presence 
of the intermediate point turned the destructive inter- 
ference into a constructive one. Changing eo to posi- 
tive (negative) values then simply shifts the whole curve 
to the left (right). Changing ei to non-zero values gen- 
erates either a "Fano" zero in the conductance (when 



eo < 0), similar to the results in Fig. 5, or a "Fano" res- 
onance (when eo > 0). At eo = one observes a change 
from a zero to a resonance as t\ changes from positive 
to negative values. All of these parameters can be easily 
changed using the various gate voltages, e.g in the setup 
of Ref. 23. Without going into much further discussion, 
it is clear that the intermediate point on the side branch 
is very effective in changing the interference pattern be- 
tween Fano-Kondo resonances and anti-resonances. Data 
from such experiments can then be used to obtain infor- 
mation on Odd- 
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FIG. 6. Conductance of the "T" network, with an extra site 
"1" between the dot and the "intersection" site "0" . Plots are 
for U = oo, T = 0, J x = .2, A b = .1 (all in units of J). From 
top to bottom: eo — 0.3, 0, — 0.3. For each value of eo, the 
figure contains graphs for ei = —0.5, 0, 0.5, with increasing 
thickness. 



VII. DISCUSSION 

We have derived an approximate analytic expression 
for a complex mesoscopic network, which contains a 
quantum dot with electron-electron interactions, and 
which may connect to several leads. Our formulae correct 
and generalize earlier expressions, and give a good quali- 
tative interpolation between the Fermi liquid behavior at 
very low temperatures and the simpler high temperature 
one. 

Although our formulae reproduce many features re- 
quired by the Fermi liquid theory, they are still approx- 
imate, and should thus only be used for discussing the 
qualitative variation of various quantities on the param- 
eters characterizing the network. However, the simplicity 
of our expressions allows for relatively easy comparisons 
with experiments and with more complicated numerical 
work. 

We have demonstrated the use of our formulae for the 
two simple cases of a single quantum dot attached to two 
leads and of the "T" network, where the dot is connected 
to another site which couples to the two leads. Indeed, we 
have reproduced and extended all the expected phenom- 
ena for these two cases. In an upcoming publication 30 
we shall apply this scheme to the interesting case of the 
Aharonov-Bohm interferometer. 
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APPENDIX A: CALCULATING T = M' 1 

We need to invert the matrix M. = u> + %r\ — H nc t- We 
start by dividing our Hilbert space into two parts, called 
"1" and "2" . Then, we write 



M 



Mn M 12 
M 2 i M 22 



and the inverse matrix T as 
T = 



F\\ Fyi 
F21 F22 



It is now easy to obtain the Dyson-like equation, 
Fn = {M 11 -M 1 2[M22]- 1 M 2i y l . 



(Al) 



(A2) 



(A3) 



In our calculations wc only need T mn for to, n within 
the finite mesoscopic part of the network (possibly in- 
cluding the dot). Identifying this part with the subspace 
"1" above, we thus only need the matrix elements of the 
finite matrix Fn. To obtain these via Eq. (A3), we need 
the Green function for the leads, [M22] -1 . For our Eq. 
(12), the matrix M22 separates for the different leads, 
and we end up with 



F?' 



M 2 



E^ 



(L)] 



with the self-energy matrices 



(A4) 



(A5) 



Here, the matrices M22 and T,^ are all of order N x N, 
where iV is the number of sites in the mesoscopic part of 
the network. Also, 



vW 



E 



(r\k)(k\s) 



(A6) 



is the Green function for the disconnected lead L. 
For the one-dimensional leads, we use the identity 



2 , JdnF n l ~ S r CT + Gdda- 
n 

— (1 + T>oGdda) 



du)' 



-E 



1 



Q *-£ cj + ii] - e k 



2iJsin \q\a 



(A7) 



with u) = — 2 J cos qa, to obtain the self-energy due to the 
lead L, which is attached at the single point til, 



£ (L) M m „ - -\JA 2 Z l]q]a 5 m .nJ m .nJJ- 



(A8) 



The same result also applies when the mesoscopic part is 
empty, and the leads are connected directly to the dot. 



APPENDIX B: DETAILS OF THE EOM FOR T (2) 

Following the same logic that led to Eqs. (22) and 
(33), the formal solution to Eq. (34) gives 

— 2_^ Jdn^nl = SoTo- — ^_^ Jdm^mn [{X n )Gdd(7 
n mn 

+ {dl a a n _ a )(l + Y> G dda )]. (Bl) 

Calculating the thermal average (d_ a a n - a ) via the 
fluctuation-dissipation Eq. (32), and using Eqs. (20) 
and (24), we end up with an integral which requires the 
sum 



5*1 — 2^ Jdm,3~{uj) m l3~{uj')i n J n d- 



(B2) 



mill 



For this and for similar sums, it is helpful to use the 
identity 



T{u + irj)f(J ± irf) 



T{J ± irf) - T{u + irf) 
u> — u) 1 + it] 



(B3) 



which follows from Eq. (26), together with a careful use 
of l/(x + irf) = V(l/x) - in5(x). Using Eq. (23), the 
sum in Eq. (B2) then becomes 



Si 



1 



LU — LU 1 + IT] 



[E (w')-SoM]- 



(B4) 



Similar manipulations allow the calculation of the ther- 
mal averages (X n ). Finally, we end up with 



/(c 



2iri ) u> — u>' + ir\ 
liti) uj — ui' + ir\ 



(l + S^ d _ ff )(S -4)-(c'^c'*) 



G'dd-ai^O - Sq) - (w' -^ U)'*) 



(B5) 



As explained after Eq. (36), Gdda and So & re understood to be functions of u> — > in- 
dependence on u>' — > u>' + ir), with u>'* — > u>' — ir\. 



it], and the primes denote a 
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